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Abstract 

The nucleon-nucleon interaction in the region of the nucleon kinetic energy up 
to 1000 MeV is analysed together with the reaction jd — > pn in the photon energy 
range = — 400 MeV. Nine nucleon-nucleon s-channel partial amplitudes are 
reconstructed in the dispersion relation N/D method: ^Sq, ^Si Di, ^Pq; ^Pi, 
^Pi, ^^2) ^^2) ^-^2 and ^^3- Correspondingly, the dispersive representation of 
partial amplitudes A^A — > pn, NN* — > pn and NNtt pn is given. Basing on 
that, we have performed parameter- free calculation of the amplitude jd pn, 
taking into account: (i) pole diagram, (ii) nucleon-nucleon final-state rescattering 
'yd ^ pn ^ pn, and {Hi) inelastic final-state rescatterings jd — > A^A(1232) — > pn, 
'jd — > A^A^*(1400) — > pn and 'yd NNtt — > pn. The jd — > pn partial amplitudes 
for nine above-mentioned channels are found. It is shown that the process ^d — > 
pn pn is significant for the waves ^So, ^Po, ^Pi, at = 50- 100 MeV, while 
"yd — > A^A pn for the waves 3P2, '^D2,^F3 dominates at E^ > 300 MeV. Meson 
exchange current contributions into the deuteron disintegration are estimated: they 
are significant at E^ = 100 - 400 MeV. 



1 Introduction 

Investigations of the lightest nuclei such as deuteron, ^He, ^H, or '^He were always chal- 
lenging for researchers. There are at least two reasons for that: 

(i) Investigations of simple nuclei give a chance to comprehend better the crucial problems 
of nuclear physics, such as the structure of nucleon-nucleon forces, the role of multinucleon 
interactions, the role of internuclear forces at small and large distances, the problem of 
multiquark bags like 6q or 12q, meson degrees of freedom in nuclei, the role of relativistic 
effects, etc. 

(m) The study of light nuclei provides an opportunity to elaborate the methods for the 
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description of composite systems. Hadron physics is a physics of composite particles — 
quarks and gluons being constituents. Elaboration of new approaches for the description 
of composite systems based on the study of light nuclei has an undoubtful advantage that 
constituents (nucleons) are observed particles, with known interaction amplitudes. 

This paper is devoted to simultaneous investigation of the two two-nucleon processes, 
namely, A^A^ ^ A^A^ at < T < 1000 MeV and -fd ^ pn at < < 400 MeV. This 
simultaneous investigation is not accidental: the A^A^ —>■ NN amplitude is basic for the 
study of the two-nucleon reactions. At the same time, the reaction jd — > pn is one of the 
simplest processes with two nucleons; therefore, the simultaneous analysis should clarify 
the points: 

(i) To what extent and at what energies the nucleon-nucleon interaction governs the other 
processes induced by this interaction; and 

(a) What characteristics of the nucleon-nucleon interaction play the leading role. 

For the deuteron photodisintegration at the photon energies in the laboratory frame 
< < 400 MeV, numerous experimental data exist, see refs. - [|13| and references 



therein. This region covers both small energies, where quantum mechanics works, and 
comparatively large energies, where relativistic effects and inelastic processes are impor- 
tant. 

The reaction 'yd — >■ pn is investigated using dispersive integration technique, which is 
appropriate for the analysis of partial wave amplitudes: 
{{) This technique is relativistic invariant; 

(m) It keeps under control all intermediate states allowing to take properly into account 
final state rescatterings which are determined by the right-hand singularities of the am- 



plitudes [0]; (Hi) In this technique there is no ambiguity related to the mass-off-shell 
amplitudes. 

The dispersion relation technique has much common with both non-relativistic quan- 
tum mechanics approach and light-cone variable technique, the methods which are widely 
used by non-perturbative QCD (or Strong-QCD Q). A similarity of the dispersion re- 
lation and light-cone variable techniques was specially outlined in ref. |0 where the 
Perturbative-QCD onset in meson form factors has been analysed. 

Consider in more detail the mechanism of the deuteron photodisintegration, 'yd — *• 
pn, at E^ < 400 MeV. The following processes contribute to this reaction: Fig. la — 
the deuteron break-up without nucleon interactions in the final state. Fig. lb — the 
deuteron break-up followed by the rescattering of outgoing nucleons. Fig. lc,d — the 
photoproduction of one or several mesons in the intermediate state, or — Fig. le,f — 
the photoexcitation of resonances, like A and A^*(1400), followed by the transition of the 
produced particles into the pn state. 

We start with the discussion of the nucleon-nucleon interaction. There is a striking 
property relevant to the hadron interaction at low and intermediate energies: in a wide 
energy range the nucleon interaction is almost elastic or quasi-elastic (this means an ex- 
citation of nucleons into highly located baryon states), i.e. with the energy growth the 
onset of genuine inelastic processes is delayed. Typical example is given by the deuteron 
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channel {pn scattering in the ^Si and ^Di waves) where, though the threshold of the 
pion production is at T = 290 MeV, the phase shift analysis indicates that the inelas- 
ticity parameter is strongly suppressed up to T = 1000 MeV |r7|. For other waves the 



inelasticities are related to baryon excitations; for instance, in the and ^F^ waves 
the inelasticities are due to the transitions NN — > A^A(1232). The same effect is seen 
in the pion-pion amplitude, where the inelasticity onset in the s-channel, / = 0, starts 
at the energy much higher than the four-pion threshold. Presently there is quite a lot 
of evidences that this is a general feature of the processes at low and intermediate en- 
ergies: inelasticities are mainly related to the transition of interacting hadron into an 
excited state, while the prompt production of new quark-antiquark pairs is suppressed. 
This property being observed several decades ago (see, for example, |]l8l and references 



therein) is discussed starting from the early 60's; it is a foundation of duality ideas and 
models, such as Veneziano's model [p!9[] . 

The direct way to perform a relativistic description of nucleon-nucleon interactions 
is to use the invariant Feynman technique that leads to the Bethe - Salpeter equation. 
However, the direct application of the Bethe - Salpeter equation to the nucleon-nucleon 
interaction faces problems related to mass-off- shell amplitudes and meson production 
processes. Mass-off- shell phenomena are usually taken into account by the introduction 
of form factors, that undoubtedly affords an ambiguity in the quantitative representation 
of the process. 

Another problem arises in the description of meson production within the Bethe- 
Salpeter equation. Indeed, the solution of this equation is represented as a sum of ladder 
diagrams of the type shown in Fig. 2. The cutting of ladder diagrams defines inelasticities 
of the amplitude; an example of such a cutting is shown in Fig. 2 by dashed line, it 
corresponds to the transition NN NN + meson. The problem is that the forces in the 
deuteron channel related to the pion exchange are rather significant. Therefore, due to the 
cuttings of Fig.2 - type, the Bethe — Salpeter equation leads to significant inelasticities 
which are switched on at invariant energies ^/s = 2m + m^^, 2m + 2m^, etc. To suppress 
the early inelasticities, in line with the suppression of genuine inelastic processes, and at 
the same time to hold a non-small interaction related to pion-exchange forces, one should 
take into account other more complicated diagrams, in order 1) to cancel the right-hand 
s-channel singularities with the help of a destructive interference, and 2) to keep the 
neighbouring t-channel singularities (or left-hand ones) that respond to the non-small 
meson-exchange forces. 

Moreover, an exact summation of ladder diagrams is hardly realizable in the realistic 
case, thus necessary simplifications being induced. One of them consists in the con- 
sideration of a potential-type interaction. Another one is an actual replacement of the 
constituent (nucleon) propagator by its residue {k^ — m^)~^ — > 2Tri9{ko)S{k^ — m^) |20[| . 
However, this latter procedure leads to anti-casual singularities of the amplitude; hence, 
one should take care these singularities be far from the region the approach pretends to 
describe. 

A correct handling of the s- and t-channel singularities of the amplitude may be 
performed comparatively easily in the framework of the dispersive N/D method 



where contributions of the left-hand singularities due to the t- and w-channel exchanges are 
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independent of the right-hand singularities which are related to the s-channel elastic and 
inelastic processes. An advantage of this technique is its relativistic invariance, together 
with an absence of a problem of the off-shell amplitude definition. For all intermediate 
particles, the equality k'^ = is fulfilled, thus allowing us to describe interactions using 
mass-on-shell amplitudes. Moreover, in the dispersion relation technique there are no 
principal problems with the description of inelasticities related both to the excited baryons 
and prompt meson production, because all s-channel intermediate states are under control. 
This means that, when using dispersive integrals, one can take into account certain set of 
the nearest singularities in due course. 

It is also important that the dispersion relation technique allows one to determine the 
relativistic wave function (or vertex function) of the composite (deuteron-like) system. 
The description of two-particle composite systems within the framework of the dispersive 



N/D method was developed in ref. [21|, where a Lorentz-covariant amplitude for the two- 
particle composite system interacting with an external vector field has been written, and 
ambiguities related to the dispersion relation subtractions have been eliminated basing 
on the gauge invariance and analyticity requirements. Therefore, while in the Feynman 
diagrammatic technique the gauge invariance is automatically fulfilled and the interaction 
amplitude obeys the Ward identity, in the dispersion relation technique these general 
properties of the theory should be considered as supplementary constraints. This is a 
necessary payment for treating independently left- and right-hand singularities. 

In ref. the dispersion relation analysis of the pn scattering amplitude in the 

^5*1 — ^ Di channels has been performed. Basing on the phase shift data, the vertices 
of the pn scattering in this channel have been constructed together with the deuteron 
vertex functions which are relativistic analogues of the deuteron wave functions. The use 
of deuteron wave function obtained in such a way allows to describe well the deuteron 
form factors up to Q"^ = 2 GeV^/c^ as well as deuteron binding energy and magnetic and 
quadrupole moments. A successful description of deuteron properties has been obtained 
taking properly into account the relativistic effects in the pra-component, without including 
non-nucleonic degrees of freedom into deuteron wave function. 

The application of the dispersion relation technique to the description of a deuteron 
has a long history (see refs. p2|): at the early stage, the deuteron form factor has been 
considered using the t-channel dispersion relation, thus bringing the t-channel anomalous 
singularity into consideration. The recent deuteron description in terms of dispersive inte- 
grals is based on the double dispersion representation in the deuteron-mass channels, 
and it does not require special treatment of anomalous singularities. In this point the 



technique of ref. |21[| is analogous to the non-relativistic wave function technique. The 
detailed consideration of the composite-particle wave functions in the dispersive technique 
and its relation to the technique of light-cone variables may be found in refs. ||2^ . 



Successful description of the elastic deuteron form factors within the dispersion relation 
technique, gives rise to the problem of the application of this technique to the other two- 
nucleon processes, in particular to the deuteron photodisintegration. 

The reaction 'jd —>■ pn was studied during several decades ||2^-||3^ in the framework of 
other techniques, mainly in the quantum mechanical and Feynman diagram approaches. 
These investigations clarify many features of the reaction, among which — an important 
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role of the meson exchange currents (typical diagrams are shown in Fig. 3). 

The analysis of the deuteron photodisintegration in the framework of the dispersion 
relation technique have been carried out in [Q-[3£] starting from small energies, where 
the contribution of inelastic processes is negligible, with a subsequent increase of the 
energy. In the present paper we expand the analysis onto the region of the photon energy 
E^i ^ 250-400 MeV, where the production of the lightest baryons, A(1232) and A^*(1400), 
is important. 

The investigation of the reaction '-yd — »• pn in the dispersive integral technique com- 
plements the investigations carried out within other approaches. This is related to the 
fact that similar processes from the point of view of their diagrammatic presentation may 
be treated in different approaches in their own ways. Such an example is provided by 
meson exchange currents, for in the dispersive technique they come from the contribution 
of anomalous t- and M-singularities to the amplitude 76? pn. 

The paper is organized as follows. 

In Chapter 2 the nucleon-nucleon amplitude is analysed in the framework of the dis- 
persion relation technique in the energy range up to T = 1 GeV. The aim of this analysis 
is to restore nucleon-nucleon amplitude vertex functions, which play the role of forces of 
the quantum mechanics approach. These vertices are used in the next Chapter for the cal- 
culation of the final state interaction of nucleons in the deuteron photodisintegration. In 
Section 2.1 the operator expansion of nucleon-nucleon amplitude in partial waves is done, 
together with the method of constructing covariant partial wave operators. In Sections 
2.2 and 2.3, the dispersion relation N/D method is described: it allows to separate left- 
and right-hand singularities or, in other words, to single out the interaction forces from 
the production processes. Here the equation for the partial amplitude is written for the 
one-channel case when inelastic processes are small, and it is also shown how to present 
the amplitude using the diagrammatic language. In Section 2.4, basing on the phase shift 
analysis, the N/D amplitude is restored in the channels ^Pi and ^-D2, where inelasticity is 
small. In Section 2.5, the diagrammatic technique is generalized for the two-channel case. 
This allows us to include inelasticity into the nucleon-nucleon amplitude and expand our 
analysis on other waves. The channels ^1)2, and with the inelasticity related to the 
production of A(1232) are considered here. Basing on the phase shift data in the channels 
A^A^ and A^A, the amplitudes are restored for the coupled channels '^D2{NN) — ^ S'2(A^A), 
^P2iNN) -5 P2iNA) and ^F^iNN) PsiNA). In Section 2.6 we briefly discuss the 
waves ^5*1 and ^Di. these waves were investigated in detail in ref. [^. In Section 2.7 
the structure of inelasticities in the non- resonance channels ^Sq, ^Pq and ^Pi is studied: 
several mechanisms are considered, such as the production of A(1232), A^*(1440) or the 
S-wave ttN pair. The technique for the calculation of three-particle loop diagrams which 
are necessary for the description of inelastic processes is discussed in Appendix A, there 
also the A^-function parameters are given. 

In Chapter 3, within the dispersion relation technique, the amplitude pn is 

calculated up to the incident photon energy = 400 MeV. The photodisintegration 
amplitude is taken as a sum of diagrams shown in Fig. 1. Dashed blocks in Fig. Ib-f cor- 
respond to the nucleon final state interaction, which includes elastic-scattering diagrams 
of Fig. 4a-type and intermediate-state particle production processes of Fig. 4b,c-type. 
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In Section 3.1 the pole diagram of Fig. la is written in terms of dispersion relation ex- 
pression in the deuteron channel. In Section 3.2 final state interaction is calculated for 
the channels ^5*0, ^5*1 Di, ^Pi, ^Pq, ^-Pi, ^-^2, ^-^2, and ^F^ taking into account the 
diagrams of Fig. 4a- type only (i.e. neglecting inelastic rescatterings) . Such an approx- 
imation allows one to describe the photodisintegration amplitude up to = 10 MeV. 
In Section 3.3 inelasticities related to the production of A(1232) are included into the 
final-state channels ^P2, and ^F^. Two mechanisms contribute to the intermediate 
A-production: the A-production in the final-state nucleon-nucleon interaction block (Fig. 
4b) and the prompt A-production by photon (Fig. le). In Section 3.4 the contribution 
of inelasticities into 70? pn is calculated for nucleon-nucleon waves ^Sq, ^Pq and ^Pi. 
It is shown that the process 'yd ^ pn ^ X ^ pn, where X is one of inelastic states 
{N7t)sN, A^*(1400)A^ or A^A, dominates in these waves. These transitions are important 
at E^ = 50 — 100 MeV, while the transition 'yd NA pn in the waves ^D2, ^P2 and 
^Pg, Fig. le, is important at E^ ~ 300 - 400 MeV. 

The detailed discussion of the results is presented in Chapter 4. 



2 Dispersion Relation Representation of Nucleon-Nucleon 
Scattering Amplitude 

2.1 Nucleon- Nucleon Vertex Operators for Partial Amplitudes 

Consider the structure of multipole operators for the NN partial scattering amplitudes. 
There are two ways for the presentation of the amplitude expanded in multipoles. Stan- 
dard matrix element for the transition A^A^ — NN (see Fig. 5a) is as follows: 

{^ip[)Q'^iPi)){^{p'2)Q'^{P2)), (2.1) 

where and are multipole operators. However, in the dispersive technique it is more 
convenient to use charge-conjugated fields: 

Vi>,(V2) = -^iP2)C , ^ci-P2) = C^{P2) , (2.2) 



where C is charge-conjugation operator. Using the Fierz transformation |^0[, "^{pi) 
\E'c(P2)? one gets 

{^{p[)Q'^,{-p',)) (^,(-p2)Q'^(Pi)) , (2.3) 

where Q'^ and Q'' is another set of operators corresponding to the "^^^^Lj NN state. 
The whole amplitude, A{s,t), depends on s = {pi + P2Y and t = (pi — p'lY and can be 
represented as a sum of partial amplitudes, A^{s), as follows: 

A{s,t) = ^(^'(p;)4.Q^...^c(-p'2))(^^c(-P2)4Qj....^(pi))A^(3) . (2.4) 

b 

Here and Qj^^,... are left and right operators of the partial amplitude, fi,^... are 

indices related to the to the state '^^'^^Lj, and I^j is isotopic operator. Summation in eq. 
(2.4) is carried out over the whole set of operators Q^^y...- 



6 



Consider the structure of the Q^^^... operator. For the spin-| particles, these operators 

are constructed using metrical tensor (yf^j,, antisymmetrical tensor e^^aB, 7-matrices and 
momenta fc^ = \{pi — P2)^i and = (pi + P2)ix- Let us introduce the tensor Tf^i^i...^ij_i,... 
related to the two-nucleon loop diagram: 

{Sp[-iPi + m)Ql,...{-p2 + m)Ql,,,.]) = (2.5) 

Here 5""^ is the Kronecker tensor, the brackets (■ ■ ■) denote the integraging over two-particle 
phase space, d$(pi,p2): 

rf^(Pi,P2) = \ ^-^^{2nmP - p, -p2)5(m'-pl)5(m'-pl), (2.6) 

where 



/d*(p,.P,) = ^^^^pW. (2.7) 

We impose the following normalization constraint: 

(2.8) 

For the states with total angular momenta J = 0, 1, 2, the T-operators are equal to: 



J = 
J= 1 
J = 2 



r = 1 



Ti.',^ = gi'^, , (2.9) 



where gf^^ = g^^ - Pf,Pu/s. 

With these requirements, we get the following set of the pn vertex operators for the 
states ^'^+^Lj with isotopic spins / = 0, 1: 



^^o(/ 


= 1) 


'Po{I 


= 1) 




= 0) 




= 0) 


'Pi{I 


= 0) 


'Pi{I 


= 1) 


'P2(I 


= 1) 




= 1) 




= 0) 




= 1) 




X 



Q 



_ 1 

i/2(s-4m2) ' 



1 

2^5 
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( 



2£(£-4m2)75^^, 
s2 (s-4m2) ^Ma/J-^^a-P/g^^ ; 



2s(s-4m2) 
n - v^60 



kj^r, + k^r,-l{k^r)g^ 



/lO 



i^v f^ti 3 9, 



HI/ 



kj^^uapSlaPph 



70s2(s-4m2)3 

Ixi/^abcxlfbPc ~l~ i/-^^ abcixlfbPc ~l~ abcflfbPc ■> Q 



Q = ^Q, 

Qfi ~ Qui 

Qij, — Qui 
Qh ~ QiJ-1 
Q/j,!/ — Qfj-i/i 
Qui/ — Q/j-i/i 
Q ^iv Q^ivi 



Q, 



Hi/X ■ 



(2.10) 



7 



liereki^k^-{Pk)P^/s, j^^^^-PP^s, = 7^ - 2A;^/(2m + v^i) and 

^n^x ~ ^M^f'^x" V(^A'^4:~'~^^^x/'~'~^xfl'/«')' stress that in (2.10) we have introduced 

left and right operators because the trace (2.5) may be the negative in the technique used: 
this leads to different left and right operators. 

It should be noted that ^'c and ^' have opposite parities, so ^db"^ is a scalar and ^c"^ 
is a pseudoscalar. 

For the singlet and triplet isotopic states, the isotopic operators are: 

/ = 0: lij-^, / = 4 = ^, (2-11) 

where Tij arc Pauli matrices. 

Consider as an example the construction of the operator Q^^, for the ^P2 state: the 
total momentum J = 2, so this operator is traceless symmetrical tensor of the second 
rank constructed of /c^ (the P wave state) and T,^ (the spin-triplet state). Therefore, 
the operator for the ^P2 state is proportional to O^i, = kj^Ty + k^T^j^ — '^{k-^V)gj^^, where 
Qnv = N{s)0^y and N[s) is normalization constant. We have: 

2 

{Sp[-{pi + m)0^y{-p2 + m)0^fyf]) = -s{s - Am )T^,.,^/^/p(s). (2.12) 
Therefore, the operator for the ^P2 state is: 



Likewise, one can construct vertex operators for all partial amplitudes. 

The partial amplitude A^{s) for a fixed state is obtained by projecting the whole 
amplitude A{s,t) given by (2.4) on this state. Corresponding conjugate operator is as 
follows: 

(*(pi)gi,...*e(-p2))(^c(-p'2)Q^.'...*(p'l))- (2.14) 

Thus, we have: 

4T^,...^,,,,„p{s)A\s)p{s) = j d^p,,p2)d^p[,p'^) 

x(^'(Pl)gt....^c(-p2))A(s,t)(^'e(-p'2)Qiv'...*(p'l))- (2.15) 

After performing azimuthal-angle integration, the following expression for A''(s) is ob- 
tained: 

d 7 — — - 

^Y(^(Pl)Q;:....^e(-P2))^(«,t)(^c(-p'2)Qtv'...*(Pi)), (2.16) 

where z is the cosine of the angle between pi and p\ in the center-of-mass system of 
particles 1 and 2: z — l + 2t/{s — Am?). 
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2.2 Left-Hand and Right-Hand Singularities and Unitarity of 
Partial Amplitudes 

In this Section we remind the analytic structure of the amphtude when nucleons interact 
by pion exchanges (m^ is pion mass). In this case the amphtude has a pole in the t-channel 
at t = m^, as well as threshold singularities at t = (nm^)^ (n = 2, 3, ...). In the s-channel 
the amplitude has a singularity at s = 4m^ related to nucleon rescattering and branching 
points s = (2m + nm.^)'^ related to inelastic channels. For the deuteron channel a bound 
state with the mass M exists: the amplitude has a pole singuliarity at s = M^. The s- 
channel partial amplitudes have the same s-channel singuliarities (right-hand ones) that 
the whole amplitude, A{s,t), has. Besides, partial amphtudes have left-hand singularities 
which correspond to the t- and w-channel singuliarities of A{s, t) with branching points 
s = 4m^ — m^'n? — see Fig. 6. Neglecting inelasticities means that we neglect all right- 
hand cuts, except for the first one with the branching point at s = Am?. 

Consider the unitarity of the partial amplitude in the physical region (the upper edge 
of the elastic cut). One has SS^ = 1 with 

5 = 1 + t{27r)^6\p, +P2- p[ - p'2)A{s, t). (2.17) 

Using the decomposition of the amplitude (2.4) and equation (2.5), we get the unitarity 
relation for the partial amplitude: 

ImA\s) = p{s)\A\s)\^ (2.18) 



2.3 One-Channel Scattering in the Diagrammatic Dispersion 
Relation Technique 

According to the N/D method |T^, the partial amplitude A{s) can be represented as a 
ratio of two functions (here and below we omit the index b in the partial A^A^ state): 

A(s) = ^, (2.19) 

with N{s) containing left-hand sinquliarities and D{s) containing right-hand ones only. 
It follows from the unitarity condition (2.18) that 

ImD{s) = -p{s)N{s). (2.20) 

The used normalization for D(s) reads: D(s) ^ 1 at s ^ oo. It is also assumed that 



there are no CDD-poles |^ in the amplitude. In this way one gets: 

0(s)^l - r 'Keim^^,.B(s). (2.21) 

^ ' JAm^ 71 S' -S ^ ' ^ ' 

One can restore the A^-function using experimental data, thus solving the problem of 
finding an amplitude with correct analytic properties in the vicinity of the point s = 4m?. 
However, to describe the nucleon-nucleon system (say, deuteron) which interacts with an 
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external field, it is convenient to introduce vertices, which, in the case of a deuteron, play 
the role of wave functions [^]. In the simplest case, the vertex function may be defined 
as follows: 



Gis) = ^Nis) . (2.22) 
Then A{s) can be written as a sum of dispersive diagrams shown in Fig. 7a. 

A{s) = G(s)(l + B{s) + B{sf + ■ ■ ■)G{s), (2.23) 

where B{s) corresponds to the one- loop diagram with the vertex G{s). 

Consider the one-loop dispersive diagram and its connection with the Feynman one- 
loop diagram. The Feynman integral reads: 

^ r d'k a{kl{P-hY-,e,{P-kf)a{e,{P-kf-,kl{P-k,f) 
^ J {2n)H {'m? - k'^){'m? - {P - kY) ' ^' ^ 

Here a{k\^ {P—kiY\ k"^, {P — kY) is an irreducible block, without two-particle intermediate 
states. The receipt for the transition from Feynman integral to the dispersive one is the 
following. We need: 

1) To neglect right-hand multiparticle singularities in the block a; 

2) To factorize the block a: 

a{kl, {P - k^Y; k\ {P - kf) G{kl {P - kif)G{k\ {P - kf)- (2.25) 



3) To calculate imaginary part of the one-loop Feynman diagram, considering the inter- 
mediate state with the energy ^/s as a real one (s > Am?): 

d^k 1 



{2t:)H {m? - k'^){rr? - {P - kf) 

d^k 1 



{2'k)H 2i 



{2mye{ko)6{m' - k')e{Po - ko)6{m' - {P - ky) = p{s); (2.26) 



4) To restore the 5-function using dispersion relation integral over the imaginary part 
related to the elastic cut, discAp: 

°° "f^') = r (2.27) 



4m2 IT S' — S Jim? IT 

Finally, we get G{s)B{s)G{s). 

The above consideration demonstrates the applicability limit for eq. (2.24): the basic 
assumption is a factorization of the block a. 

A general method to perform a factorization of the t- and w-channel meson exchange 
interactions has been developed in [Q. For the description of a composite system in 
general case, one should introduce several vertex functions Gi related to a separable 
interaction Vf. Ni ^ Vi = GiG^, where the left vertex function, Gi, and the right one, G\ 
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may be different. For further consideration, it is convenient to introduce the amphtudes Oj 
which do not contain the right vertices G". The whole amplitude and auxiliary amplitudes 
tti are related as: 

A = J2aiG". (2.28) 



The amplitudes Oj obey a set of linear equations: 



where (see Fig. 7b) 



ds'&{s')p{s')Gi{s') 



s' — s 



Rewriting (2.29) in a matrix form, one has 

d = Bd + g , 

where 

d = 

The final expression for the partial amplitude reads: 

A = f{I-B)-^g, 

where 



ai 
a2 


, 9 = 


' G' ' 


, B = 


\B\ 
Bl 


Bl .1 

Bl . 















(2.29) 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



r = (2.34) 

The amplitude obtained in such a way is unitary; it may be derived from the unitarity 
condition directly. 

Now let us explain how the scattering amplitude discussed above may be written in 
terms of the energy-off-shell Bethe-Salpeter equation. We need to introduce the energy- 
off-shell scattering amplitude: 



A{s,s^) = Y,a^{s)G\s^). 

i 

Using equation (2.31), one can write this amplitude as 

ds' p{s' 



(2.35) 



A{s,si) 



A{s,s')V{s', si) + y(s,si). 



(2.36) 



Um? 71 S' — S 

This is a dispersive Bethe-Salpeter equation for ^4(5, si), where the effective interaction is 

V{s,si) = J2Gi{s)G'{s,). (2.37) 

i 

The suggested diagrammatic technique is a modification of the standard N/ D method: 
this technique is suitable for the description of an interaction with an external field. 



11 



2.4 One-Channel Approach for the Waves ^Pi and 

For the waves ^Pi and the inelasticity is small, and corresponding amplitudes can 
be described in the one-channel approximation. We restore vertex functions of the pn 
scattering, relying upon the phase shift analysis data |T^. The G-function is determined 
by the left-hand singularities, therefore it is represented as a dispersive integral along the 
left-hand cuts: 

Gis) = r (2.38) 

J-oo TT S — S 

where sl = 4m^ — is the location of the nearest pion branching point. Below, to 
simplify a cumbersome calculation, we replace the integartion in (2.38) by the summation: 

G^.{s) = j:^^, si = 4-h\j-l), (2.39) 

with the poles placed on the left from sl- 

The values of parameters 7*, and Sq obtained by fitting the phase shift data are 
discoursed in Appendix A. The results of the phase shift fit together with the experimental 
data are demonstrated in Fig. 8. 

2.5 Two-Channel Scattering, NN and NA: 
the Waves ^Da, and ^F^ 

In this Section we consider the inelasticity in the nucleon-nucleon amplitude due to 
the production of the isobar A(1232). The nucleon-nucleon amplitude is treated as a 
two-channel one (the channels A^A^ and NA). We analyse here the coupled channels: 

(A^ AT) -5 ^2(A^A), ^P2iNN) -5 P2(A^A) and ^F3(A^A^) PsiNA). 

Matrix element for the A^A^ —>■ NA transition (see Fig. 5b for the momentum notation) 
is written in the form 

(^"(p'i)Q"A,^o^c(-p'2))(^c(-P2)g^;,/^(pi))A"^(5). (2.40) 

Here and are vertex operators for the states a and b. T\ and I are isotopic operators 
in the A^A and A^A^ channels, correspondingly. Aab{s) is the partial amplitude for the 
transition stateb state a. \E'^ is the wave function for the spin-| particle; \E'^ obeys the 
equations: 

p^^^ip) = 0, 7m^m(p) = 0. (2.41) 



For the A^A states, vertex operators follows: 



'S2{NA) : Q 



^P2{NA) : Qal3,fi — £abcakaPbSc(3,fj. + £abcl3kaPbSca,fi 
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2 



(2.42) 



Here fli is a normalizing factor; the momenta are defined as P = pi + p2, k = — P2), 
s — P^, /c* = — P2) — ^(^^A ~ rn'^)P, where uia is the A-isobar mass; A;^, 7^ and 
g^p are introduced in Section 2.1, while is defined as 



2fc^ 



m + niA + y/s 



(2.43) 



The A^A^ vertex operators. Q^. for ^D2, ^P2 and states are given by equation (2.10). 
The isotopic operator Tq for the states A+n and A°p is defined as 



1 

71 





1 
1 




V2 



10 
10 



(2.44) 



The partial amplitude >1 is 2 x 2 matrix. It is related to the 2x2 ^'-matrix as follows: 



S = I + 2i^A^, 
where y/p is the matrix of square roots of phase space factors: 





-/PnX 



As before (see Section 2.3), the partial amplitude is represented in the form: 

A = f{I-B)-'g, 
but with another structure of vertex operators: 



Gjv Gf 
Ga 



(2.45) 



(2.46) 



(2.47) 





\ 


" 




G* 





9 = 





G" 







G^ 



(2.48) 



and S-matrix 





Bnnn 





BNNt 


" 


B = 


BtNN 





Btm 













BrAA 







BaAv 





Baaa . 



Bijk{s) 



ds'G'{s')pj{s')Gk{s') 



s — s' 



(2.49) 
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Here G^G^, GaG^, G,.G* and GtC are A^-functions for the transitions NN NN, 
NA A^A, NN NA and A^A NN, correspondingly (see Fig. 5e), Ti = Am'^, 
T2 = (2m + m^)^. In our fit we put G^ = Gr and G* = Gt] as it was in previous Section, 
tlie two- vertex form is used for the transition A^A^ A^A^: GatG^ = GiG^ - G2G2 — 
Fig. 5f. 

For vertex functions the parametrization (2.39) is used. Parameters for G are found 
using nucleon-nucleon and nucleon-isobar phase shifts, 6nn and 6n/\, and the inelasticity 
parameter t] which are defined as 



Su = r/e^*^-- ,\S,2\ = \S2l\ = ^Jl-v' , S22 = ve''' 

Here indices 1 and 2 refer to the channels A^A^ and A^A. The description of the phase 
shift data in the coupled channels ^Z?2 S2, ^P2 P2 and — ^ P3 is shown in Fig. 9 



"JVA 



together with the experimental data [43 . 



2.6 Coupled channels ^5*1 and 



In Fig. 8 the coupled channels '^S'l — Di are drawn, for which the inelasticity is small, 
though the connection of channels is important and the two-channel approach must be 



applied here [El 



2.7 Models for the Inelasticity in 5*0 wave 

According to Jl^, the mechanism of inelasticity in the waves ^Sq, ^Pq and ^Pi differs from 
the inelasticity in the channels ^D2, ^P2 and ^Fs, and the amplitude has a non-resonance 
behaviour. 

Here we consider several models for inelasticity in the ^5*0 wave. First, we suggest that 
the inelasticity comes due to the production of A(1232). Second, we consider as a source 
of inelasticity the production of A'^J'^(1440). This resonance is heavier than A(1232), but 
it has larger width, and the transition amplitude ^Sq{NN) — ^ Sq{NN*) is not suppressed 
near the threshold. Third, the prompt pion production in the ttN S-wave is assumed. 
It should be noted that the 1st and 2nd variants correspond to the production of the 
TrA^-state in the P-wave . 

1) A(1232)-isobar production 

Matrix element for the transition ^So{NN) — ^ Do{NA) is determined by (2.40). The 
vertex operator for the ^Dq wave is 

Kai3 ,Q^ = -Qm' (2-50) 
where Q is the normalization factor and 

K^^ = kik'^-lk''9^^. (2.51) 

The other quantities refer to equations (2.42) and (2.43). 

2) Production of A^*(1440) 



'Do{NA) : 
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Matrix element for the transition NN — > NN* is written in the form similar to (2.40): 

(f^^*(p;)g«rvi/,(-p'2))(^c(-P2)g'r^(Pi))A«'(s). (2.52) 

The resonance A^*(1440) has the same quantum numbers as a nucleon, so the vertex 
operator for the NN* state is: 

^So{NN*): g = ^75 = (2.53) 

V 



3) Production of the (7rA^)5-pair in the J = ^ state 

Since the lightest resonances in the -kN system, with quantum numbers = \ , are 
located far from the pion-nucleon threshold, we deal here with the case of the prompt 
non- resonance pion production (Fig. 5c). Nevertheless, we approximate the ttN pair by a 
remote quasi-resonance R, with an energy-dependent width. This enables us to operate 
with two-particle intermediate states using the same technique as was used for A (1232) 
and iV*(1440). 

Consider the vertex of the quasi-resonance decay R —>■ ttN: = c^r. The width of 
the quasi-resonance R, F^, is defined by its decay into pion and nucleon. This allows us 
to relate the quantities c,r and F^. The resonance propagator is written as: 

(p + Mr) 



Ml-p^-B{p^)' 



(2.54) 



where B{p^) corresponds to the pion-nucleon loop diagram (Fig. 5d). Real part of the 
loop diagram re-determines the mass of the quasi-resonance Mr, and imaginary part of 
B{p^) provides its width. Let us rewrite the propagator (2.54) as a series: 

MI-p^-IMrTr {Ml-p^)^ {Ml-p^f ^(M|-p2)^- ^-"'^ 

where the first term corresponds to the propagator of a stable particle and the second 
one to the propagator with the inclusion of the pion-nucleon loop diagram. The Feynman 
integral for the second term of the right-hand side of eq. (2.55) is equal to: 

(j^^Mh) f d'k 4(/'- + "0 (j> + Mb) ^2.56) 



Ml-p"^ J {2t:)H {w? - k^){ml - {p - kf) Ml-p^ ' 

The comparison of eqs. (2.55) and (2.56) provides the expression for F^ in terms of c^. 
To this end, let us decompose k in external vectors p and k-^, the latter being orthogonal 
to p: 

k = ap + k^, (2.57) 

where a = (p^ + m^ — m^) /2p^. The term which is proportinal to k^ vanishes in eq. (2.56). 
When comparing eqs. (2.55) and (2.56), it is necessary to calculate imaginary part of the 
loop diagram and to neglect its real part; for this purpose the following replacement should 
be done: 

(m^-k^)(rll-{p-kr ^ - ky{ko)5{ml - {p - kmpo - k,) . (2.58) 
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Final expression for the width takes a form: 



R 



SvrM, 



R 



[s + m 



m. 



2\2 



.1/2 



4s 



:2.59) 



Vertex operator for the NR system in the ^5*0 state is: 

1 



^SoiNR) : Q 



1, Q = Q. 



(2.60) 



The description of the A^A^ amphtude is performed within the framework of the two- 
channel amplitude, where the first channel is the A^A^-state and the second is A^A, or 
NN*, or NNtt states, depending on the model we use. The amplitude is restored basing 
on the phase shift analysis data of ref. [O. The results of fitting to the phase shift 



analysis data for the ^5*0 wave are shown in Fig. 10a,b. All the three models work equally 
well in the description of the experimental data on the elastic amplitude A^A^ — > A^A^. 
However, they provide different results for the amplitude A^A^ — NNir. 



2.8 iViV*(1400) Channel as a Source of Inelasticity in ^Pq and ^Pi 
Waves 

The description of the ^Pq and ^Pi NN waves is performed assuming the production 
of A^*(1400) as a source of inelasticity in the intermediate state. The results of fitting 
procedure for the ^Pq and ^Pi waves are shown in Fig. lOc-f. Fitting parameters are given 
in Appendix A. 

3 Deuteron Photodisintegration within the Disper- 
sion Relation Technique 

In this Section the deuteron photodisintegration amplitude is calculated basing on the 
analysis performed in the previous Section. This amplitude is represented as a sum of 
the pole (impulse approximation) diagram (Fig. la), diagram with nucleon-nucleon final 
state interaction (Fig. lb) and diagrams with the photoproduction of A(1232), A^*(1400) 
and pions in the intermediate state (Fig. Ic-f). It should be stressed that all parameters 
which are necessary for the calculation of these diagrams have been found in the Section 
2, where nucleon - nucleon interaction have been studied; therefore the calculation of the 
■yd — > pn amplitude is parameter- free. 



3.1 Pole Diagram 

Consider the amplitude A^°'', which corresponds to the pole diagram of Fig. 11a, written 
within the dispersion relation technique. We start with the presentation of the Feynman 
pole diagram: 

AP°' = ^^e,MrJ, (3.1) 
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where 

M-' = *(Pi)r,|^/^r^*e(-P2) . (3.2) 

Here and are polarization vectors of the deuteron and photon which are orthogonal 
to the deuteron and photon, P (P^ — M^) and q [q^ = 0), momenta: 

(eP) = 0, {eq) = Q. (3.3) 

The deuteron isotopic operator, 7'^(7 = 0), is 

4- I- (3.4) 

Vertex functions for the transition deuteron — > nucleons, F^, and for the photon-nucleon 
interaction, F^, are of the form: 

Ff = 7.0+(p-p2).V^, = 6(7^/1 + + Pi U2), (3.5) 
where the vertex functions /i and /2 are matrices in the isotopic space: 

fi-iff + f-n + iff-fnrs. (3.6) 

Here 1 and T3 are isotopic matrices 1 = diag{l, 1) and T3 = diag{l, —1). Indices p and 
n stand for the photon-proton and photon-neutron form factors; 0, </? and fl^{N — p,n) 
depend on invariant variables. 

Consider the transition from the Feynman pole amplitude (3.2) to the dispersion 
relation one. For this diagram the dispersion relation is written in the deuteron channel, 
where the pole diagram contains the two-particle threshold only, which corresponds to the 
production of two nucleons in the intermediate state (in Fig. 11a the relevant two-particle 
cut is shown by dashed line). 

When going from the Feynman amplitude to the dispersive one, the intermediate 
state vector P should be introduced, see Fig. 11a: P = p + p2, P^ = s. The transition 
to the dispersive integral in the deuteron channel is carried out at fixed total energy 
* — (Pi +^2)^- Let us introduce q — —p + pi which satisfies the equations: 

P + q = Pi+P2 = P + q- (3.7) 

For the calculation of kinematic relations, it is suitable to use the deuteron rest frame, 
where 

q={qo:q): P=(M,0), P=(VI,0), q = {M + q^ - ^/§,q). (3.8) 

One can see that q^ — {p — Pi)^ 7^ q^- This inequality is related to the energy non- 
conservation in the dispersive integrals. The photon energy in the deuteron rest frame go 
is determined by s as: 

+ 2Mqo = s. (3.9) 

We get from (3.9) and (3.10): 



f^s+~s-j^{s + M''). (3.10) 
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Calculation of the discontinuity in the deuteron channel is equivalent to the replace- 
ment of the propagator by the ^-function, considering the two-nucleon intermediate state 
as a real one: 

^ h2m)Sim'-p')9ipo). (3.11) 



Using the discontinuity across the two-nucleon cut, we reconstruct the amplitude 

= / , ^^^(^' - p')0{po)^{pi)T,{p + m)I''Ttmci-P2) . (3.12) 

The substitution (3.11) corresponds to physical decay of the virtual deuteron with momen- 
tum P into two nucleons with the on-shell momenta p and p2- Rewriting the (5-function 
argument in the form 

p"^ - = s - 2pp2, (3.13) 
we elimitate the integration in (3.12): 

= ■^^^{Pi)'^f^{P -P2 + m)I^ri^,{-p2) . (3.14) 



s 



M2 



In the photon-deuteron center-of-mass frame one has: 

^"V^ AT') '- = ^°<^'^0' 

Let us compare the expressions for the pole diagram provided by the Feynman ampli- 
tude and the dispersion relation one. In the deuteron rest frame the denominator of the 
Feynman amplitude (3.2) is: 

2(m2 - = 2(2Mp2o - M^) = 2M(2p2o - M). (3.16) 

Corresponding denominator related to the pole term in the dispersion relation amplitude 
is equal to: 

5 - = Aplf, -M^ = (2p2o + M)(2p2o - M). (3.17) 

Let us stress that the momentum p is different in both techniques, namely, p = P — p2 
for the Feynman integral and p = P — p2 for the dispersive one and, correspondingly, 
the functions 0, (p, fi and /2 are different too. However, in the non-relativistic limit, 
(2p2o + M) ~ 2M, the dispersive and Feynman pole amplitudes coincide with each other. 

Since the intermediate particles in the dispersive technique are mass-on-shell, the 
functions (j) and in (3.5) depend on P^ only (P^ = s). Their numerical values are given 
in pri. Form factors and depend on only. 

These functions can be expressed in a standard way in terms of the electric, , and 
magnetic, nucleon form factors: 



In equation (3.8) is equal to zero. 
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3.2 Nucleon-Nucleon Final State Interaction 



We calculate here the triangle diagram of Fig. lib, with nucleon-nucleon interactions 
in the final state. In the triangle diagram the nucleonic degrees of freedom are taken 
into account only; contribution of non-nucleonic degrees of freedom is considered in the 
next Section. Although we consider here the one-channel nucleon-nucleon amplitude, the 
developed method can be easily generalized for the case of multichannel scattering in the 
final state, such as the scattering in ^5*1 — ^ Di waves. 

We classify partial 7(i pn amplitudes according to the final state pn-waves {pn 
system in the state b is determined by (2.4) and (2.33)). The amplitude 'yd — > pn with 
the pn final state interaction is written as 

Al = ^^e,M'^^, (3.19) 

where 

M',^ = f{l - 5)-1A.^,„,^,(^'(pi)4Q^;3^c(-P2)) . (3.20) 

The amplitude A^-^q,/? corresponds to the triangle diagram contribution. Fig. lib, while 
M^^ is a sum of graphs with a full set of the final state nucleon rescattering diagrams. Ver- 
tex matrix and loop diagram matrix B are determined by the pn scattering amplitude, 
see eqs. (2.30) - (2.34). A^^ ^/^/ can be written in the form analogous to g: 



V '■■ J 



(3.21) 



As in the calculation of the pole diagram, we start with the Feynman amplitude for 



^ I d% -Sp{Til\p' + m)T,{p\ + m)Ga^ph{-p'^ + m)) 



After summing over isotopic indices, one gets 

Sp{l'T,h) = iff + (-l)Vr)7. + (/I + {-iyf2)iP + P2), ^ r,- (3.23) 

To make a transition from the Feynman amplitude (3.22) to the dispersion relation 
one, it is necessary to carry out the following steps: 

1. To expand the numerator of (3.22) in external four-vectors, for which we choose P 
and Pf = pi +P2: 

-Sp {Til\p' + m)f,(p\ + m)G,Qi^l\-p', + m)) = ^ d)Fl^^^^ . (3.24) 

i 

Here F^^^^^ is an expression built up from the vectors P,Pf and tensors gTfj.,£^iTa/3- The 
procedure for constructing invariant functions rf* and the explicit form of ^.^ for various 
pn scattering channels is given in 
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2. Using (3.24), to represent the expression for ^\-^^af3 ^^e form: 

K,,.p = llD]Fl^,.p. (3.25) 

3 

For the presentation of D'- as dispersion relation amphtudes, the double discontinuity of 
the triangle diagram is calculated, i.e. the intermediate two-particle state is considered 
as a real one, with the neutron-proton energy a/I before and \fs' after the interaction 
(s, s' > 4m^). This procedure corresponds to the replacement: 

^ {27Ttf6{m^ - p'^ip'o) (3.26) 



(m^ — p'2)(m2 — Pi){7ii? — p'2) 4 

x5(m2 - p^)e{p\Q)6{m^ - P2)0{p'2o) = p{s, s, z) . 

3. To reconstruct the function using the dispersion integral with a double disconti- 
nuity as an integrand, : 

■ /""^ ds ds' disc disc 
" 7r2 (5-M2)(s'-s) ■ 

Here 

disc discD^j = J ^2 • (3.28) 

Using the results obtained by the derivation of (3.14), we perform integration over s in 
(3.27). This gives: 

0-=f^rU^),J3^M^- (3.29) 



TT 7-1 2 ' (s - M2)(s' - s 

p(s') is given by (2.7), s is defined according to (3.15), s = Pj, z = cos{q,p\). Invariant 
functions determine the function A,-^ (3.25) and, consequently, the amplitude Mf^. 

3.3 Photodisintegration Amplitude with A(1232) in the Inter- 
mediate State 

In this Section we include into the photodisintegration amplitude the A-isobar produced 
in the intermediate state. The importance of an isobar in the deuteron photodisintegration 
was emphasized in numerous papers (see, for example, ||33| , , |^ ) . The isobar reveals 
itself in two ways: 

1) As an A^A intermediate state in the final state nucleon-nucleon rescatterings (dashed 
block in Fig. lb,c); 

2) In the prompt photoproduction, with a subsequent transition NA — * pn in the triangle 
diagrams of Fig. le-type. 

Here we calculate the photodisintegration amplitude with A^A(1232) in the interme- 
diate state, taking into account the coupled channels ^D2{NN) —^S2{NA), ^P2{NN) 
-^P2{NA), and ^FsiNN) -^Ps^NA). 



20 



1) A-isobar, partial width and the coupling 

We start the treatment of A(1232) in the deuteron photoproduction amplitude with 
a consideration of the transition coupling for A^7 — »• A, which determines the process of 
Fig. le. The propagator of the stable A-isobar, with mass mA, is equal to 



m\ — 



(3.30) 



Re-definition of the propagator for the case of a decay into nucleon and photon consits in 
a replacement: 

{ml - p'r' ^ {ml - / - B,{p')r' , (3.31) 

where B^{p'^) is the photon-nucleon loop diagram. The real part of this diagram re- 
determines rriA, and the imaginary part of B^{p'^) gives the partial width for the decay 
A ^7Ar: 



(j3 + mA)(-c/^^ + |7a7^) 
ml — p^ — zmAF^(p2) 



(3.32) 



Thus, we have to calculate B^{p). The following expression stands for the A-isobar 
propagator with B^{p) considered as a perturbative correction: 

J- 7. TTT J- r- 



m. 



p. 



{27i)H 



m 



2 _ J.2 01^' _q2 



X 



(p + mA)(-gi-^/ + 
ml — p^ 



(3.33) 



The vertex operator F^^ has the form ^ 



= C^l5{q9fia - qal,^)V2To . (3.34) 
Let us decompose the vectors k and q over vectors p and k-^ {k-^ is orthogonal to p): 

q = Ap + k^ , k = Bp-k^, (3.35) 

where 



A 



1 1 
p — 



B 



p? + 



2p2 ' 2p2 
To calculate the imaginary part of (3.33) the following replacement should be made: 

- \{2^^mm- - ky{kM-qy{qo) . (3.36) 

When integrating over /c-*-, the terms proportional to fc"*" and k-^^ vanish, while the 
quadratic terms should be replaced as follows: 



k-^k-^ 



(p2 — m^) 
4p2 



9al3 — 2 ^±9al3 ■ 



(3.37) 
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Using the equality 



9ae ~^ 2 7« 7^ 



one obtains: 



m\ — 



X I 



■p ( A'^Bp'^ 



m 



-Akl 
3 ^ 



-Bkl 
3 ^ 



(3.38) 



Fixing = in the numerator, we have for the width F^: 



-pw 

^ ~ Svr 



m. 



2m/ 



m/ 



-'7 • 



After presenting as = {e/mT,)C, we get C = 0.30 (at F^^/FJ^^ = 0.006 and F 



tot 



(3.39) 
= 115 



MeV), in accordance with refs. |0 and 

2) pn-production in the waves ^D2{NN), ^P2{NN) and ^F^iNN) 
The amphtude of the process 'yd —>■ pn is given by eq. (3.20). For the waves ^D2{NN), 
^P2{NN) and ^F^i^NN)), the structure of the triangle diagram At-^q,^ is the same as g 
(which is given by (2.48)), namely: 



(3.40) 



We have two types of triangle diagrams: with nucleons in the intermediate state, A^ 
and A* (Fig. lib), and with A^A intermediate state. A'' and A'^ (Fig. 11c). The calcu- 
lation of A* is given in Section 3.2. The scheme of calculation of A* is presented below. 
The Feynman integral corresponding to the diagram of Fig. 11c is: 



A^ 





A* 













A^ 



d% -Sp{Vii\p' + m)V%{p\ + m^)P,,.G,Qlp^^,U-p'2 + m)) 
{2ti)H {m^ - p'^){ml - pf){m^ - p'i) ' I - J 

where F^^, is the vertex for the A- photoproduction (3.34). The spin-dependent numerator 
of the isobar propagator P^^' is: 



Pes' = -gi' + ' 



(3.42) 



where 



9ee' - 



P'lePls' 
Pi 



Is 



Is 



P'lP'le 
P? 
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The expression (3.41) is written for a stable isobar. In the case of an unstable isobar 
decaying into pion and nucleon, the pole term should be replaced as: 



where the function B^^lpf) is determined by the pion-nucleon loop diagram. The real 
part of this diagram re-determines the mass, and the imaginary part gives the width 
of the decay A Nit: 

Using the Lehman representation, the isobar propagator can be written as follows: 



{ml-p['-tm^T:{p['))-'= I 

J u 

where the weight function, R{fp?\ has the form: 



(m+m„)2 — p'^^ — iO 



p^-2^ 1 rnAr-(m^) 

^ ' 7r(mi-m2)2 + (^^r-(m2))2 • 



Below we accept = F];^^, and the parametrization of ref. is used for the 

momentum-dependent width: 

Y-^(m^\ - F- ^V(l + ^') 

*°*xi/(i + xi)' 

Here X = qr/m.,r, Xa = Q'A'^/'^Tr with r=0.081, q and qa are pion momenta in the isobar 
center-of-mass frame for the isobar masses rfi and m^, correspondingly. 

Transition from the Feynman amplitude (3.41) to the dispersion relation one makes it 
necessary to follow the prescription given in Section 3.2. Then we get for Dy. 

~ 9 ^/ ~ 9x f°° ds' dz 2d] 



D) = r dm^R{m') H — ( 

J(m+m„)^ J(m+m)'^ TT J- 



(m+m^)2 J(m+m)'2 TT 7-1 2 {s — ]VP){s' — S f) 



1 /2 

(s' — (m + m)2)(s' — (m — m)2)\ 1 

/ IGvr 



X ( ^ ^-^^ ^ ^ ) — . (3.44) 



3.4 Deuteron Photodisintegration and Inelasticities in the Final 
State Waves ^^o, ^Pq and ^Pi 

In Section 2.7-2.8 the description of the NN scattering in the waves ^5*0, ^Po and ^Pi has 
been done, and three different models for the mechanism of inelasticity being suggested: 
production of A(1232), A^*(1440), or (7rA^)5-pair. It is shown in Appendix A that the 
transition amplitude NN NN does not depend on the model used for the inelasticity, 
i.e. on the type of the transition amplitude NN NR, NR NN and NR NR, 
where R = A, N*,or {7iN)s. 
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This means that the process of Fig. lb does not depend on the inelasticity type as 
well, and it may be calculated within the technique presented in Section 2. However, the 
processes of Fig.l c,e,f, with a prompt photoproduction of a resonance or pion - nucleon 
pair, may depend on the inelasticity type. For the calculation of processes of Fig. lc,e,f, 
it is necessary to know the resonance photoproduction amplitudes. It should be noted 
that A^*(1400) has an about 50% branching rate into the channel Ntttt, so the triangle 
diagram with A^*(1400) in the intermediate state effectively describes the process of Fig. 
Id. 

1) Triangle diagram with A^A intermediate state 

The vertex A has been presented in Section 3.3. The calculation of triangle 

diagram with A^A intermediate state is carried out in the same way as for the waves 
considered in previous Section. 

2) Triangle diagram with A^A^*(1400) intermediate state 

The resonance A^*(1440) has the same quantum numbers as a nucleon, = | , and 
the vertex function for the transition A^7 — N* is equal to: 

rf = c^l,r, , (3.45) 

with the following relation between c-y and the width F^'of the decay A^* 'jN: 



2(Mjv. - m)2 



r^;^^* = ' ::\' ^ {Ml, - m") . (3.46) 



The correlation between and c-y may be found in the same way as it has been done in 
previous Section for A(1232). To define c^, we rely upon the fact that the total width of 
A^* is equal to Ttot = 200 MeV. Let the ratio FI^/F^^ be the same for p7 and ny channels, 
being equal to 0.10% [0. Using (3.46), we get Cj = 0.3. Experimental data on A^* 
do not allow to define the coupling with a high accuracy. 

3) Triangle diagram with N(tcN)s intermediate state 

Following the model presented in Section 2.7, we approximate a pion-nucleon pair 
by the quasi-resonance R. Then the vertex function for the transition 7A^ R can be 
written in the form: 

r? = <7M75r3 . (3.47) 
The coupling is found by calculating the 7p n'^n cross section: 

/■I 1 1 

ct(7p TT+n) = / dz—— [{s^p - {ml + m^) f - 4m^m^]^/^ 

J-i AJ Ibirs^p 

2 2 

X (M|-.,^P + M|r^^ ^^^'^^^^^ + mhMk + k + Mn) (3.48) 

2 
3 

where q and k are photon and neutron momenta, s-yp is the 7p energy squared, and 



X (A; + m)7^/75(A; + K + Mr)}- 



J = 2{s-m), g^^, = g^^, - — ^ , = g^, - k^-^ 
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The parameters c^r, Mr and Tr standing for the quasi- resonance are introduced in Sec- 



tion 2.7. According to ^8[, the 7p 7r~^n cross section has a resonance pick related to 
the A-isobar and a smooth background. We relate this background to the quasi reso- 
nance production of the (7rA^)5'-pair. The quasi-resonance width is calculated with the 
simplest ansatz about the vertex of the decay R — > nN: 

T^ = c^T. (3.49) 

The fit provides the following parameters = 0.5, c-,, = 0.7, Mr = 1.85 GeV. 

For all considered waves (^5*0, ^Pq and ^Pi) and for all types of the models, the result 
is the following: processes with inelastic production in the triangle-graph intermediate 
state provide small contribution, in other words, the contribution of diagrams Ic, le and 
If is not large, while the main one comes from the diagram of Fig. lb. For example, the 
maximum value of the triangle graph contribution in the wave ^5*0 comes at = 30 — 50 
MeV, and it does not exceed 5% of cri^So) (see below). 



4 Results and Discussion 

4.1 Polarization Averaging and Cross Section Calculations 

Within the used normalization, the differential unpolarized cross section has the form: 



\A 



2 



da = ^d^, (4.1) 

where A is the whole disintegration amplitude, which is built up of the amplitude A = 
j^poie _|_ Iqj averaging over the spins of incident particles and summed over spins of 

final particles. The summation J2b is performed over all considered states ^5*0, ^5*1 — ^ Di, 
3Po, 'Pi, 'P2, and ^F^. 

The two-particle phase space factor, d^, and the flux factor J are defined as: 

d<l> = ^-^^{27rrS\P -p,- pMm' - pl)6im' - pj) , (4.2) 

J = 4{{PqY - m^y/\ 
For non-polarized deuteron we need to make the replacement: 

^r^r' - -igi, = -I [grr' " ^) • (4.3) 



For the non-polarized photon, it is necessary to average over the two polarization states: 

-2^MM' - ~2 Y^^' P2 ^ ) ' y^-^) 

where = q — P{Pq)/P'^, and q^gj;^ = and g^J~g^^ = g^^- It means that the 
I-'--'- A , 



combination g^i^Aui which is used for the calculation of the cross section is gauge invariant. 
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4.2 Photodisintegration Cross Section: Nucleonic Degrees of 
Freedom 



At the first stage of our investigations ||39|, we liave calculated the amplitude 'jd pn 
taking into account nucleonic degrees of freedom only: the photo disintagr at ion amplitude 
is described by the sum of diagrams la and lb and parameters of G functions are restored 
on phases of A^A^ scattering only. 

Fig. 12a presents the calculation result obtained in such approximation: solid line 
stands for the pole diagram contribution (Fig. la), while dashed line describes the cross 
section a{'~fd pn) calculated with both diagrams la and lb. 

One can see that for nucleonic degrees of freedom the photodisintegration reaction is 
described at i?^ < 10 MeV only, that agrees with the result of ref. [37]. 

One should note that the diagram of final state interaction does not change the result 
obtained in the framework of the impulse approximation. This is a consequence of the 
completness condition for two sets of wave functions: plain waves and those with the pn 
interaction. The proximity of solid and dashed lines is a criterion of self-consistency of 
our calculations. It should be noted that in the region of very small the situation 
changes: the diagrams of Fig. la and lb cancel each other, due to the orthogonality of 
the deuteron wave function and the wave function of continuous spectrum. 



4.3 Contribution of Inelastic Channels 

Inelastic processes are important in the waves ^P2, and ^F^ (the production of A- 
isobar and in the waves 

^^0, ^^0 and ^Pi (the production of A^*(1400) or (7rA^)s-state). 

The calculation result for a{'~fd —>■ pn), with the account of inelasticities in the waves 
^Ps, and ^Fg, is shown in Fig. 12b (dashed-dotted line). At ~ 300 MeV the 
calculated cross section has a bump, but it is lower that the experimental one. At 10 < 
E^ < 100 MeV the dashed-dotted line is located below the curve which is given by the 
pole diagram, that is due to the destructive interference. It should be noted that in these 
waves the main contribution comes from the diagram of the prompt A-isobar production 
(Fig. le), i.e. from the process 7A^ — > A^A. 

Inelasticities in the waves ^5*0, ^Po and ^Pi lead to the considerable increase of the cross 
section a {•yd — >• pn) 10 < E^ < 200 MeV: because of that, the cross section acaicild — > pn) 
satisfactorily describes the experimental data up to E^ ~ 60 MeV (solid curve in Fig. 12b). 

Fig. 13 demonstrates the calculation results and the experimental data [10-13] for 
differential cross sections at E^ = 20, 60, 90 MeV. An agreement is seen at E^ = 20, 60 
MeV. However in the region E^ > 60 MeV there is no agreement with experimental data. 
One may believe that this discrepancy is due to meson currents which are not completely 
taken into account in our calculation. 

Let us discuss in more detail the structure of processes which are taken and are not 
taken into account in the approach under discussion. In terms of the analytic structure 
of the amplitude 'yd — >• pn, the diagrams with rescatterings do take into account the (s, t) 
and (s,m) anomalous (Landau) singularities (remind that t = {q — pi)"^ and u = (<? — ^2)^)- 
Figs. 14a,b give an example of diagrams with such singularities: Fig. 14a demonstrates 
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the diagram with the Landau singularity in (s, t)-sector, while the diagram of Fig. 14b 
has the Landau singularity in {u, t)-sector. Both these diagrams have also singularities of 
triangle diagrams, which are shown in Fig. 14c, d — these diagrams come to being due to 
the shrinkage of the upper nucleon line into a point. The diagrams of Fig.l4c-type have 
anomalous singularity of the logarithm type at t ~ m^ + 2/i^ + 4/iA/me = ttr, that is rather 
close to the physical region of the reaction (e is deuteron binding energy). Similarly, the 
singularity at m = ttr is inherent to the diagram of Fig. 14d. With final state interaction 
accounted for in our approach, we also take these singularities into consideration, but not 
completely, for just the same diagrams are present in the {u, t)-hox diagrams — see Fig. 
14e,f. These latter are in fact the meson exchange current diagrams which are not included 
into our approach, and one may believe that the lack in the cross section description at 

> 50 MeV should be filled in by the account of these diagrams, together with the 
subsequent final state interaction (Fig. 14g). 

The role of meson exchange currents in the deuteron photodisintegration has been 



estimated by Laget |3T|]. Fig. 15 demonstrates the differnce Aa = a^xpild pn) 



(^caicild — > pn) obtained in our calculations and the evaluation of meson exchange current 
contribution given in ref. for the processes shown in Fig. 14e,f. One can see that 
the estimation given in ref. 0] provides a reasonable description for Act in the region 
= 50 — 200 MeV but cannot explain the data at Ej ~ 300 MeV. This is exactly the 
region where the bump in the cross section corresponds to the A isobar in the intermediate 
state (see Fig. Ic). Therefore, it is natural to suggest that in the considered case we meet 
analogous situation: the A isobar is produced in the process induced by meson exchange 
currents — see Fig. 14h. The next transition AA^ — NN (Fig.l4i) would provide a 
contribution affecting the enhancement of the cross section around E^ ~ 300 MeV. We 
believe that this is the only reasonable explanation of the discrepancy observed in this 
region: the bump in the cross section is undoubtedly due to the A^A channel, and the only 
one which is not accounted for in our approach is the contribution of the [u, t) singularity 
into the amplitude jd — > A^A, i.e. the contribution of meson currents into the process 
-yd NA. 



5 Conclusion 

We have carried out the dispersion relation analysis of the nucleon - nucleon scattering 
and deuteron photodisintegration amplitudes. As a result, all the waves which are not 
small in the region Tiab < 1 GeV are written within the despersion N/D representation. 
A^-functions which are the analogues of the potential are restored for the waves ^5*0, 
^5*1 — ^-Di, ^Pq, ^Pi, ^Pi, ^P2, ^D2, and ^Fs (the A^-functions for the deuteron channel 
^Si and ^Di were found previously in |]2T|). Besides, the N/D representation is written for 
inelastic amplitudes A^A^ A^A (the waves ^D2iNN) ^2(A^A), ^P2iNN) P2(A^A) 
and ^F3{NN) PaiNA)) and A^A^ ^ N{NTT)s-^ave (the waves ^S'o, ^Pq and ^Pi), thus 
giving a complete dispersion relation representation of the low and intermediate energy 
amplitude at T < 1 GeV. 

The amplitude of the photodisintegration jd — > pn has been analysed in the framework 
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of the dispersion relation technique for < 400 MeV. The spectator mechanism of 
the photodisintegration related to the interaction of photon with one of the deuteron 
nucleons is considered, the inelasticity related to the production of A, N* or {ttN)s 
pair being included. Furthermore, using the constructed nucleon-nucleon amplitudes, 
the final state interaction is properly taken into account within the dispersion relation 
technique. All these calculations are parameter- free, the necessary characteristics of the 
nucleon-nucleon amplitude are directly taken from from the experiment, that is, 
amplitudes, deuteron vertices found in [^] from the fit of the deuteron form factors at 
< 1.5 (GeV/ c)^ and the NN amplitudes found in this paper. Using the language of the 
dispersion relation representation, meson exchange currents are the contributions of the 
(m, t)-singular part of the amplitude Aijd pn). The calculations performed allow one to 
find out the meson exchange current contribution: the difference between the observed and 
calculated cross sections is a direct measure of it. The analysis shows that the spectator 
mechanism for the photon-deuteron interaction, with a consistent account of the final state 
interaction, adequately describes the data at < 50 MeV. However, at E^ > 50 MeV 
the meson exchange currents begin to dominate, providing more than a half of the cross 
section at E^ > 100 MeV. Comparing the results of our calculations with those of Laget 



31[| , one may conclude that the stardard mechanism for the meson current contribution, 
namely, interaction of incident photon with a meson which provides an exchange between 
deuteron nucleons, may explain the value of the deuteron photodisintegration cross section 
at E^ < 150 MeV only, but not around 300 MeV. It favours the idea that the bump in 
Ctotild — > pn) at E'-f ~ 300 MeV with the width about 300 MeV is due to contribution 
of the (m, t)-singulare part of the amplitude 'jd — >• NA: that is a meson exchange current 
contributing into A^A^ NA when the photon interacts with t- or w-channel meson. 

At the same time, our analysis gives rise to the problem of meson exchange currents 
in the deuteron channel, the A^A^ waves ^5*1 and ^Di: the matter is that deuteron form 
factors v4(Q^) and B{Q^) are well described in a broad interval, < < 1.5 (GeV/c)^, 
without an exchange current ||21|. It is possible that in this respect the deuteron channel 
is a singular one, and meson exchange currents are not large there. It is a well-known fact 
that in the deuteron channel the inelasticities are strongly suppressed: we cannot exclude 
that in the deuteron channel these characteristic features are somehow related to each 
other. However this problem needs especial investigation that is beyond the scope of the 
present study. 

Acknowledgement We thank V.V.Anisovich, L.G.Dakhno, V.A.Nikonov and A.V.Sarantsev 
for numerous discussions and useful remarks. 



6 Appendix A 

Here the G-function parameters are displayed, which were obtained by fitting phases and 
inelasticities of the considered waves. The phase volumes are also presented which were 
used in the fitting procedure. For the deuteron waves ^Si and ^Di the parameters are 
given in ref. P^ . 

The simplest way to find the G-functions is seen by eq. (2.22). The A^-functions which 
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contain left-hand singularities related to the ^-channel meson exchange have the meaning 

of a potential of the free-nuclcon interaction. Because of that, they may change sign. This 
fact has been taken into consideration by introducing two functions Gi and G2, which are 
determined as follows: 

N = GiG^ - GsG^ (A.l) 
In the case when inelasticity appears, the A^-function takes a matrix form, namely, 

/ Nnn--^nn Nnn^na j (^2) 

the change of sign is allowed for A^tvat^atat only. 

By fitting the phase shift data, the following circumstance have been taken into con- 
sideration. In eqs. (2.45) and (2.49) the S- and i?-matrices depend on the product of 
functions G and p: GiPjG^. For the A^A^-scattering without inelasticity, pj is defined by 



eq. (2.7) pj = l/167ry'(s — 4m2)/s = po- But for practical use the phase volume factor 
Pj may be taken in such a way that, on the one hand, it would reflect the kinematics of a 
chosen wave and, on the other, it would not violate the convergence of the integral (2.49). 
Let us introduce auxilliary functions p^""^^ and cf^^^^ which are related to pj and Gj as 
follows: 

Then Gi — Gj-""^^(p^""^^/po)^^^- During the calculations the functions Q^^^^^ (not Gi) 
were parametrized as follows (2.39): (^(""^^ = Y^^=i -zrr- The introduction of is 
especially suitable for the waves with inelasticities, when one needs to include the inelastic 
threshold into the phase space. 

Now consider the partial waves. For the waves ^Pi and we use the following 
expressions for 

: ^ J^i^ (A4) 

Parameters of (^(""^^ are given in Table 1. 

The next step is to include resonances [ttN) into consideration. The modified way 
of finding out parameters is as follows. First, let the resonance width be zero: = 0, 
where R is one of A, A^* or (jrN)s- Then the phase space factor, for instance for isobar 
Pq, related to the loop shown in Fig. 4b, is equal to 



Po=J rf^AAT = — y . (A6) 

The trace over the loop in Fig.4b with the partial operators O'^^^^ in the vertices is equal: 
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Here O^^^p are partial operators which are defined by formulae (2.42), (2.50) and (2.53), 
without normalizing factor l/y/U^. P'^'^' is given in (3.42). 

Consider now the real resonance, with ^ 0, decaying into meson and nucleon, and 
instead of the diagram of Fig. 4b we use that of Fig. 4d. There arc three mass-on-shell 
particles in the intermediate state, that leads to the integral over the resonance mass for 
the phase volume, together with the normalization factor Pa{s) written in (A. 7). Then 
the phase volume for the partial state a has the form: 



and 

R(§) = ^^^^ (A9) 

""^'^ (5-M2)2 + (r-M«)2- ^^-^^ 

Normahzing factors in (2.42), (2.50) and (2.53) are equal 

= PA,a 

Now let us get the equation (A. 8): the three-particle phase space factor related to the 
diagram of Fig. 4d. It has a form: 

X (27r)^5(P -ki-k2- k^) j d%2S^{ki2 - ki - k2) J dm^d {P-k^-k2- k^) . {A.IQ) 
After multiplying it by a unity, 

j d^ki25'^{ki2 - ki- k2) j dmi25{ki2 - mi2), (^-H) 
we get a product of two two-particle phase volumes: 

They are equal to (mi = m^r, m2 = mjv, ^77.3 = mjv): 

1 dQ^N 

K-ttN' 



87^^/m?2 
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2(27r)3 v"i^ W (27r) 



1 (i^^AfA 

=A;jvA- 



(s - (mAT + mi2)2)(s - (mAf - mi2)^) rfl^ArA 



Stta/s V 4s Att 

Then, one has for the three-particle phase space factor: 



(^.13) 



I dmi2 1 dn^Adn^N /A-,A\ 



The phase space for the loop diagram of Fig. 4d is 

Here is given by (A. 7). Then 

_ r dmu Kn kNA d^lNA dVt^N / . 

87r^/;;^87rViK2-mi)2 + (rXmA)2^" 47r 47r " ^ " ^ 



(mAT+mx)^ V 12 

The parameter a may be found, using boundary condition ^ ^ata/IStt-v^) at F*" — > 0, 
providing 

c = (A17) 

^ (m^ - (mAT + m^y){m\ - {uin - m^f) 

Coming back to eq. (A. 3), one can find G = G""^(p2r''/PA)^/^ for any wave a. However, 
in this case the triangle diagram A^^^^ (3-41) should be calculated with non-normalized 
operators O^^^, for the normalizing factor has been already included into (p^T^/pa)^^^- 
Comparing cq. (A. 8) with (3.43) and (3.44) for A^^^^, one may see that the Lehman 
expansion takes explicitly into account the resonance widths, when calculating the graphs 
lc,e,f. 

In our calculations the following parametrization of the phase space factors were used 
for the waves ^D2, and ^F^: 



167r \ s s"* 



1 ^J{s - Am^y 
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2 



P2 ■■ pT = TTT- OA^ ^ — + 0.3 



IGtt \ 



1 Y'(s-4m2)3 



s — (2mjv + rriT^Y) nr^ ~ (2mjv + 2m.j^^'^' 



F3 : p^"^ = — 0.35^ ^ V ^ + 0-95 



IGtt V 



1 ^(s - 4m2)' 



For the waves ^Sq, "^Pq and ^Pi the phase volumes (p^"^) are related to eq. (A. 8) as 
pwux^g^ _ p^(^s)/s^, and i?(s) is chosen as follows: 

^W-g tS.^yVtrrm.)^ . ™.= 4.05 +(.-1)0.1. (A21) 

The meaning of such a definition of R is connected with an unclear inelasticity nature 
in the waves ^5*0, ^Pq and ^Pi. When fitting phase space data, we used as a model for 
inelasticity the set of 10 resonances with masses rrii and large widths Vf, these widths 
defined by the decay vertices to nucleon and meson, c-j^r. Ctt^ = 0.028 + {i — 1)0.028, and 



-16^'-^'^^'^'^'^''' -"^-Ji ■ ^^-^^^ 

The derivation of this formula is similar to that of Section 2.3. Parameters of G^""^) are 
presented in Table 1. 

In the ^5*0 wave there were considered several inelasticity variants. To get vertices Gi 
which respond to A(1232), N* or (7rN)s, one should calculate correct phase space pA, Pn* 
or P{NTr)s using the equations (A. 7) and (A. 8), with the mass Mr and width related 
to considered variant. Then G-functions are calculated using eq.(A.3). 
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Figure 1: Diagrams responsible for the deuteron disintegration process. 




Figure 2: An example of the Behte-Salpeter ladder diagram with nucleon rescatterings: 
the cutting of the diagram presents the connected inelastic processes. 
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Figure 3: Deuteron photodisintegration diagrams with meson exchange currents. 
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Figure 4: Loop diagrams which correspond to the elastic rescattering (a), intermediate- 
state production of the A isobar (b), prompt pion production (c), loop diagram used in 
the calculation of the two-particle phase space, with the account of the isobar width. 



37 



(2m+^)' 



Figure 6: Complex s-plane and the position of the partial amplitude singularities: (from 
right to left) right-hand side singularities at s = 4m^, (2m-|-m^)^ the position of pole, 
and left-hand ones at s = Arm? — m^, Aw? — Arn^ ... 



A(s) = X + >0< + >00<+ 








P-k, P-k P-k; 



Figure 7: Dispersion relation loop diagrams: (a) representation of the partial amplitude 
as a set of loop diagrams, (b) one- loop diagram with right {Gi) and left {G^) vertices: 
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Figure 8: Results of the fit for the phase shifts, S^n, of the nucleon-nucleon scattering. 
Tiab is the kinetic energy of incident proton in the laboratory system. The waves: a) ^Pi, 
b) ^D2] c) coupled channels ^Si —^Di. phase shifts in the waves ^Si and ^Di, Snn, mixing 
parameter e. 
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Figure 9: Results of the fit for phases and inelasticities in the two-channel approach. 
Coupled channels ^D2{NN) S2{NA), ^P2iNN) P2iNA), ^F^iNN) PsiNA): 
phase shifts in the waves: a) ^D2{NN), Snn] b) ^S2{NA),5na'i c) the transition amplitude 
squared ^ D2{N N) S2{N A) , (l-jy^); d) ^P2{NN), e) ^P2{NA), f) the amplitude squared 
^P2{NN) -5 P2{NA); g) ^Fs{NN), h) 'P^iNA), i) ^Fs{NN) Ps{NA). 
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Figure 10: Results of the fit for the phases and inelasticities in the two-channel approach 
for the waves ^5*0, ^Pq, ^Pi. a) phase shifts in the wave ^5*0, Snn, b) inelasticity in the 
wave ^5*0, (1 — ?7^); c) phase shifts in ^Pq, d) inelasticity in ^Pq; e) phase shifts in ^Pi, 
f) inelasticity in ^Pi. 
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Figure 11: 



43 
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Figure 12: Total cross section, Ut, for the deuteron photodisintegration: a) the con- 
tribution of the impulse approximation diagram la to the at (dashed line), at with FSI 
(nucleon degrees of freedom) (solid line); b) the contribution of the impulse approxima- 
tion diagram (dashed line); at with the inelasticity in the waves '^D2{NN), '^P2{NN), 
''^F3{NN) (dot-dashed line); at with the inelasticity in the waves ^5*0 (A^A^), ^Po{NN), 
^Pi{NN) (two dot-dashed hne); the complete at (solid line). is the photon energy in 
laboratory system. 
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Figure 13: Differential cross sections ^(6) at = 20 MeV, = 60 MeV, E^ = 
95 MeV. 
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